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The high positron production rate required to explain the flux of 511 keV gamma rays from the 
galactic center has inspired many models in which dark matter creates positrons. These models 
include the annihilation of light dark matter and scattering of dark matter with excited states (ex- 
citing dark matter). We show that existing cosmic microwave background data robustly constrains 
such models when the annihilation or scattering cross section is not velocity suppressed depending 
on the model of the galactic dark matter halo. Upcoming data from the Planck mission can exclude 
the fiducial Via Lactea II halo model, which also provides a good fit to the 511 keV morphology. 
We additionally find combined constraints on exciting dark matter scattering and annihilation and 
update constraints on the lifetimes of dark matter excited states. Finally, we apply constraints to 
models of dark matter decay in which produced positrons fall into the galactic center and produce 
the 511 keV signal on their annihilation, demonstrating that most of the parameter space of interest 
is ruled out. 
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I. INTRODUCTION 

A narrow line of 511 keV gamma rays from the galactic 
center has been observed since the 1970s [l|, 0], indicat- 
ing annihilation of > 10^'^ electron-positron (e='=) pairs 
per second. This emission has been studied extensively 
by INTEGRAL/SPI since 2002 The spectrum in- 

dicates that ^ 97% of the e"*" annihilate through positro- 
nium formation, which limits the e"*" injection energy to a 
few MeV or less, and a more detailed analysis finds that 
a significant majority of the annihilation occurs in warm 
neutral and ionized phases of the ISM. Morphologically, 
the emission has two components, from the galactic cen- 
ter bulge and from the galactic disk. Depending on the 
model used, the bulge/disk ratio for annihilation rates 
is - 1.5-6. 

Over the years, many models for the astrophysical pro- 
duction of mildly relativistic e"*" have been proposed, al- 
though standard mechanisms have difficulty reproducing 
the observed morphology of the emission, particularly the 
large bulge/disk ratio; [ll| gives a thorough review of 



not only the emission itself but also many e"*" produc- 
tion mechanisms. One which has received a great deal 
of attention is decay of radionuclides produced in 
supernovae and heavy stars, which reproduces the disk 
emission well. [12', 13] have argued that, in a particu- 
lar model of e+ transport, radionuclide decay can also 
source the bulge emission. On the other hand, as the re- 
view [ll| notes, the model of [H, [l^l relies on a number 
of apparently arbitrary assumptions about the interstel- 
lar medium and rnagnetic field (and its turbulence) in the 
bulge. Recently, [lj| found that low-energy positrons an- 
nihilate close to their production sites in a range of prop- 
agation models, so, in contradiction to 0, [l3l, they find 
that radionuclides (produced mainly in the disk) cannot 
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source the bulge emission. In any event, it is clearly pre- 
mature to claim that standard astrophysical mechanisms 
can produce the galactic bulge positrons. 

Due to the lack of consensus on candidate astrophysical 
sources for the galactic bulge emission, more exotic mech- 
anisms have been of great interest in the past decade, in- 
cluding a number involving dark matter (DM) . Roughly, 
DM models for the 511 keV excess may produce pairs 
through the decay or scattering of DM particles. There 
are two cases of DM decay that produce e'^ with only 
mildly relativistic injection energies: first, the DM itself 
can have an MeV-scale mass [H, [l^, or the DM can 
have several states with the deca y th rough an MeV-scale 
mass gap producing an pair [17f. In either case, as- 
suming small propagation after production, the 511 
keV flux is proportional to the integral of pdm along 
the line of sight. As the DM density is expected to in- 
crease toward the galactic center, all these models nat- 
urally have a large bulge/disk ratio for production. 
However, [H, ^ showed that the INTEGRAL signal is 
more highly peaked toward the center of the galaxy than 
realistic DM density profiles (even the most cuspy found 
in simulations) , essentially ruling out decay of DM at ei- 
ther mass range as the source of the galactic bulge e"*" 
excess. Recently, f20l] proposed a novel model in which 
a small fraction of DM decays, producing e"*" outside the 
galaxy; these e"*" fall into the galaxy and annihilate (pref- 
erentially in the galactic center, as claimed in (20l]). 

Alternately, DM can produce pairs through scatter- 
ing processes. In the simplest such models, the DM has 
an MeV-scale mass, and the scattering process is direct 
annihilation to an pair. These models and additional 
signatures have been studied extensively in [l^ [2l| - [3l| . 
Alternately, more massive DM with several states and 
MeV-scale mass gaps can scatter into an unstable excited 
state, which decays into the ground state via emission of 
an pair, as discussed in 17, 32-43]. (These are excit- 
ing dark matter, or XDM, models.) In these scattering 
models, the gamma ray flux follows the line of sight in- 
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tegral of so reasonable halo models can produce 

the observed signal; in fact, [T^ has found that the halo 
model derived in the Via Lactea II simulation [13, 
has a maximum likelihood ratio very close to the peak 
value. 

The lack of a strong astrophysical candidate for the 
galactic bulge emission as well as the surprising (and 
striking) agreement between the emission morphology 
and simulated DM halos motivates us to search for other 
potential signals of DM scattering models either as cir- 
cumstantial evidence in favor of these models or as con- 
straints on them in the case of non-observation. For ex- 
ample, the XDM models discussed in [H [H [H, EO, El] 
naturally include a several hundred MeV gauge boson 
with weak coupling to electric charge, which could be 
discovered at fixed target experiments such as the Mainz 
Microtron [46] or APEX In this paper, we place con- 
straints on XDM and light annihilating DM explanations 
based on limits from the cosmic microwave background 
(CMB) anisotropy spectrum; as we will review, energy in- 
jected into the Standard Model (SM) plasma around and 
after the era of recombination is tightly constrained by 
large I CMB anisotropics. As long as the XDM scattering 
or light DM annihilation cross section is not suppressed 
at low velocities, we find that these models for the 511 
keV emission will be constrained by forthcoming results 
from the Planck satellite [i^, specifically ruling out the 
preferred DM halo parameters from Via Lactea II. 

The plan of this paper is as follows. In the next sec- 
tion, we systematically review DM models for the pro- 
duction of mildly relativistic in the galactic center. 
Then we present a brief review of existing constraints on 
DM annihilation (or scattering) from CMB anisotropics. 
In section ITVl we present constraints on creation dur- 
ing (and after) recombination and discuss how they con- 
strain XDM and annihilating MeV-mass DM models. 
Next, since XDM models contain decaying excited states 
of DM, we discuss how CMB observations provide con- 
straints on those lifetimes, updating the results of 
and also place constraints on the recent model of 
We conclude with a brief discussion of our results. 



II. DARK MATTER MODELS FOR THE 511 
KEV EMISSION 

As discussed above, the morphology of the 511 keV 
emission rules out models of DM decay as an explana- 
tion for positron production in the galactic center, so we 
will focus on models in which DM scattering processes 
produce pairs. The rate of positron production in 
such models is 
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{(7Vrel){x) , (1) 



where M is the DM mass, s is 1/2 for real/Majorana 
DM and 1/4 for complex/Dirac DM, rj is the number 
of e+ produced per scattering event, and Y = Y/Ydm 



is the relative abundance of the active DM state, a 
is the appropriate scattering cross section with rela- 
tive velocity v^ei, and (• • • ) is the average over the DM 
velocity distribution at position x. We further write 
(avrei) — o^{F{vrei)) , whcrc CTw is the cross section at 
some representative relative velocity and -F is a dimen- 
sionless model-dependent function of the relative veloc- 
ity. The e+ production rate is i? = 1.1 x 10^^ s'^ 
We follow [13] in writing 



R - Anr^sY^pl — ( (kpc)^ , 



(2) 



where Pq is the DM density in the solar neighborhood, 
which we take to be Pq = 0.4 GeV/cm^. ( is 



C = kpc-M drr'ipDMir)/pQy{F{vrei)) , (3) 
Jo 

where we take the bulge radius Vc — 1.5 kpc, correspond- 
ing to the width of the INTEGRAL signal. We assume 
that the DM density follows an Einasto profile 

PDMir)/pQ = exp [~{2/a) {{r/r^r ' (^o/^s)")] (4) 

with the sun located at =8.5 kpc and a and as free 
parameters. The Via Lactea II simulation is fit by a = 
0.17, rg = 25.7 kpc, which also lies in the best fit region 
of [l^ for the 511 keV signal. Smaller values of a and 
lead to a more cuspy central halo. We assume that DM 
velocities follow a Maxwell distribution with dispersion 
Wo(r) cut off at the escape velocity Wesc('') satisfying 

t;o(r)3 (X r^-^VW , wo(r'0) = 220 - 230km/s (5) 



v,sc[rf = 2v^{rf [2.39 + ln(10 kpc/r)] 



(6) 



as in [l3] (the choice of escape velocity follows [49], and 
the velocity dispersion is s ugg ested by simulations includ- 
ing baryonic contraction [50|). 

We will consider three models in which DM scatter- 
ing processes produce mildly relativistic positrons. The 
first is annihilation of MeV-scale DM, and the following 
two are endothermic and exothermic XDM respectively. 
Finally, we will review a recent model by f20| in which 
decays of a metastable DM component creates e^ at late 
time which fall into the galactic bulge. 



A. Annihilation of Light Dark Matter 

As first proposed by [2l|, DM particles of mass 2ra^ < 
^5 10 MeV have a sufficient e+ production rate if 
they annihilate into pairs with the appropriate cross 
section, av ~ 10~^^(M/MeV)^ cnf/s. As this cross sec- 
tion is several orders of magnitude smaller than required 
for the correct thermal relic abundance, either the an- 
nihilation is p-wave dominated or annihilation to is 
subdominant to another channel, as in [2^. Morpholog- 
ical studies of the 511 keV signal using a similar profile 
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Contours of log (ann (dashed) and log (^_|_ (solid) 


as 



function of the Einasto profile parameters a and Vs. The dot 
represents the Via Lactea II parameters, d is calculated for 
vt = 10"^ 

for the velocity dispersior0 disfavors p-wave annihilation 
[T^ , so we consider only the s-wave cross section. In this 
case, F(vrei) = 1 and C, is simply the integral of p%]^j- 
Figure [1] shows contours of log C (henceforth denoted as 
Cann for this modcl) as dashed lines as a function of a 
and r-g. 

A number of authors have considered additional pre- 
dictions of this model. [1^ [H, [2^ studied the predicted 
511 keV emission from nearby dwarf spheroidal galaxies, 
which is consistent with observed limits (these predic- 
tions should be similar for XDM models). [1^ demon- 
strated that models like those of 28] with additional light 
states could be detected in neutrino oscillation experi- 
ments. Recently, [13, HH demonstrated that DM with 
MeV masses annihilating preferentially into as op- 
posed to V can alter Big Bang nucleosynthesis, providing 
an alternate set of constraints. Finally, [1^ [23| placed 
constraints on MeV DM annihilation to based on the 
CMB; our results extend and update those works. 

B. Endothermic Exciting Dark Matter 

XDM as first envisioned by [H, IH, HH] consist of a DM 
ground state (state #1) and an unstable excited state 
(state #2) with a mass splitting SM12 > 2me- DM colli- 
sions above a threshold velocity vt populate the excited 



state, which then decays to the ground state by releas- 
ing an pair. With Af ^ TeV and SM12 > 2me, vt 
is approximately the velocity dispersion of DM in our 
galaxy. In order to reduce vt or to accommodate lower 
DM masses (to ^ 5 — 10 GeV), consider that a signifi- 
cant fraction of DM remains in a metastable excited state 
(state #3), and upscattering of that state into the unsta- 
ble state through a smaller mass gap SM32, allowing a 
sufhcient e+ production rate [13, [la, [13, l39l-|4ll. Direct 
detection of XDM has been considered in |l7l. l37l - [39j . 

In these models, the fraction Y of DM in the 
metastable excited state is determined by freeze out of 
DM-DM scattering, including the same scattering re- 
sponsible for present-day e+ production, and is further 
influenced by the kinetic coupling between DM and SM. 
Details for models with gauge kinetic mixing are provided 
in [17]. One consequence of XDM is that a potentially 
significant fraction of DM is initially (after chemical and 
kinetic freeze out) in the unstable excited state. [13] has 
provided limits on the lifetime of this decay based on 
CMB measurements; our results in section |V] will update 
those constraints. 

In these models, F{vrei) = \/'v1eih't - 1 0(«re! - Wt), 
including both a phase space suppression and the kine- 
matic threshold. Since Vt ^ 10""^, the kinetic energy of 
DM at recombination and later (before structure forma- 
tion and virialization) is much too small to allow this 
upscattering process to proceed. Therefore, endothermic 
DM models are kinetically forbidden from producing 
pairs during the era of interest and are not subject to the 
CMB constraints discussed in section HVl below. 



C. Exothermic Exciting Dark Matter 

Another possibility is that the metastable DM excited 
state #3 scatters exothermically into the unstable ex- 
cited state #2 [l7|; in fact, this possibility is preferred 
when the DM is charged under an abelian gauge group 
[4^ . This downscattering process is not kinematically 
suppressed, so it can produce p airs through out recom- 
bination. In this case, F{vrei) — V^^re/Z^t' + l[3 Very re- 
cently, [43 proposed an exothermic XDM model that also 
produces a gamma ray line at the DM mass, simultane- 
ously explaining a line at 130-135 GeV in galactic center 
observations of the Fermi satellite [sH - lssf (see [s^ for a 
review, including some references to DM models for this 
signal) . 

In exothermic XDM, we define the velocity Vt as the 
threshold velocity for the inverse upscattering process (or 
equivalently, as the velocity imparted to the less mas- 
sive states for downscattering at rest). This is given by 
ILvf/2 = KSM23, where fj. = M/2 is the reduced mass 



Our velocity dispersion decreases slightly less rapidly as r — > 
for a fixed density profile. 



^ There is a weak additional velocity dependence; details for scat- 
tering by gauge boson exchange are given in llTll . 
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and K is the number of DM particles that excite or de- 
excite in the scattering (typically k = rj). For example, 
for non-Abelian XDM as discussed in [131, k = 2 and 
vt = 2y/26M23/ME For example, a mass of M = 10 
GeV and vt = 10~^ correspond to a mass splitting of 
1.25 keV between the two excited states of DM. The solid 
contours in figure [T] show contours of log C (henceforth Cj^ 
for exothermic XDM) as a function of the Einasto pa- 
rameters a and for Vt — 10^'^. Compared to Cann 
a fixed a,rs, Ci is enhanced, and the enhancement in- 
creases as Vt decreases. As a result, contours of log^i 
shift to the right as vt decreases and approach the con- 
tours of log (ann as Vt oo. Sincc Vt < 10"'^ typically 
requires additional fine-tuning, we consider the contours 
of figure [1] to give an estimate of the uncertainty in our 
modeling. In addition, constraints placed on models with 
Vt — 10~^ will be weaker than models with larger vt, so 
we consider it to be a conservative choice. 

The other phenomenology of exothermic XDM is sim- 
ilar to that of endothermic XDM. In some cases, the 
metastable state can decay to the unstable state by emis- 
sion of a single X-ray photon; this signal can be near 
observable levels flTII. 



D. Infalling Positrons from Decaying Dark Matter 

Recently, "2^ proposed a new model of decaying DM 
which could provide a sufficient number of e+ to explain 
the 511 keV signal. In their model, DM has two states, a 
stable ground state and a metastable excited state with 
mass splitting 6M < GeV. The excited state can decay 
to the ground state by emitting e+, and the constraint 
on the mass splitting prevents antiproton production (al- 
though a more general model may be acceptable from 
that point of view). The e"*" cool to nonrelativistic ener- 
gies via scattering processes; some in large orbits are 
only entering the galaxy and annihilating in the present 
day. Assuming the ground and excited states are similar 
in mass, [20i] found that Yrec/M ~ 5 x 10"* GeV"^ yields 
a sufficiently strong gamma ray signal (as above, Y^ec is 
the fraction of DM in the excited state at recombination). 
This model is free from the morphological constraints of 
[isl [l9| because the DM decays occur outside the galaxy, 
but no detailed morphological study has yet been per- 
formed. In section [Vj we will find constraints on this 
class of models independent of details of the decay. Due 
to the presence of an excited DM state, we will include 
these models in our definition of XDM as a form of short- 
hand. 



III. DARK MATTER AND THE COSMIC 
MICROWAVE BACKGROUND 

Any energy injected into the SM inter-galactic medium 
(IGM) from a hidden sector after matter-radiation equal- 
ity modifies the recombination history of the universe in 
two ways. First, the injected energy can thermalize and 
heat the IGM, or it can directly ionize hydrogen and he- 
lium atoms. The overall effect of the additional energy 
is not to delay recombination but to increase residual 
ionization, increasing the optical depth of the IGM and 
thickness of the last scattering surface, which slightly al- 
ters both the TT (mostly suppressing the high multipole 
moments) and polarization anisotropy spectrum of the 
CMB. A standard assumption of the weakly interacting 
massive particle (WIMP) paradigm for DM is that the 
DM annihilates into SM particles with a cross section set 
by the thermal relic abundance. As a result, there is a 
wide literature deriving constraints on excess energy in- 
jection based on CMB measurements, particularly due to 
DM decay or annihilation. For example, see [1^, [13, [stI - 

The physical quantity that determines any modifi- 
cation of the CMB anisotropy spectrum is dE/dVdt, 
the rate of energy deposition to the SM IGM per vol- 
ume. It is common to separate astrophysical factors 
from model-dependent particle physics; specifically, we 
can parametrize the standard cases of DM decay com- 
pletely to SM states 



dE 



dVdt 



DM decay 

and annihilation into SM particles as 



dE 



dVdt 



plnl^ii + zfpa 



(7) 



(8) 



DM annihilation 



Here, pdec — f /t is the DM decay rate multiplied by an 
energy deposition efficiency /, discussed further below; 
note that, if all DM decays, the lifetime must be long 
compared to the age of the universe. Also, 



Pa 



2fs- 



M 



(9) 



Note a factor of 2 redefinition of vt compared to |17|| . 



where Uann is the total annihilation cross section and 
other variables are defined as in section [ll] above. Con- 
straints on the modification of the CMB anisotropics 
from standard cosmology translate into constraints on 
Pdec and Pann ■ These constraints currently exclude some 
annihilation channels for thermal WIMPs of mass < 10 
GeV; they can be weakened in asymmetric DM models 
[66l mf\ or avoided if the DM dominantly annihilates into 
hidden sector particles (68j . 

In the case of XDM models, we are also interested in 
decay of an excited DM state to the ground state plus 
e^e~ . In this case, the pair carries off a fraction 
5M12/M of the DM mass energy. Accounting for the 
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fact that only the unstable excited DM states decay, we 
find pdec — if /t){SMi2/M)Yu, where the fractional 
abundance of DM in the unstable state, also depends 
on redshift for lifetimes r shorter than the age of the 
universe due to the exponential decrease of the decaying 
population. Exothermic XDM scattering is similar to the 
annihilation case in that it involves 2 DM particles, so we 
replace Pann in (|H]) with 



Pscatt = Y^fsr]6Mi2 



M2 



(10) 



with additional variables defined in section III CI Note 
that endothermic XDM scattering is kinematically for- 
bidden at recombination. Constraints on Pann translate 
directly into constraints on Pscatt- 

The deposition efficiency / depends on the specific par- 
ticle physics model of DM through the branching ratios to 
different decay/annihilation products and their eventual 
cascade into charged SM particles, photons, and neutri- 
nos. At high energies, / is well approximated by 1 minus 
the fraction of energy carried off by neutrinos. However, 
in addition, since the IGM is transparent to photons of 
certain energies at certain redshifts, / depends on red- 
shift z, as well. Several articles have studied the redshift 
dependence of the deposition history in detail, notably 
including , and [TOI presented a principal component 
analysis delineating how the redshift dependence of pdec 
and Pann determine the CMB anisotropics. This redshift 
dependence has been included in the newest CMB con- 
straints on decaying and annihilating DM as reported in 
[(35, 67, 71- 74]. In particular, J2] noted that the full red- 
shift dependence is well approximated (within 15%) for 
DM annihilation by taking / to be constant at the value 
for z — 600 (which is also the center of the first principal 
component of [t^). For moderately relativistic e^, [g^l 
finds f{z = 600) > 0.9, and we adopt that constant value 
(the error associated with the approximation is smaller 
than the allowed range in 5M for XDM models). 

We will compare energy deposition from XDM mod- 
els to updated constraints from the CMB, specifically 
constraints that are independent of annihilation chan- 
nel. In particular, |71i] places a 2(T constraint of Pann < 
5.73 X 10-27 cm3/s/GeV using WMAP7 data [zl and 
forecasts a constraint oipann < 9.78 x 10~^^ cm^/s/GeV 
from Planck. Similarly, [T^] finds a 95% confidence con- 
straint of Pann < 2.42 X lO'^^ cm^/s/GcV from WMAP7, 
Pann < 2.09 X 10-27 cm3/s/GeV from WMAP7 and ACT 
[Tq . and a forecast constraint of Pann < 3.03 x 10-^^ 
cm^/s/GeV from Planck. Finally, 'ts'] find pann < 
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7.86 X 10-28 cm^/s/GeV using WMAP7 and SPT 
data at 95% confidence. We will also discuss the decay 
of XDM excited states and the decaying DM models de- 
scribed in section III Dl While we reserve the details of 
constraints on decaying DM for section |V] below, we note 
that [13, [S^] agree on a 95% confidence limit of roughly 
Pdec < 2.4 X 10-2^ s"^ for lifetimes longer than the age 
of the universe. 

Finally, we note that energy deposition in the SM IGM 




FIG. 2. Contours of log Pann (in cui /s/GeV) for s-wave an- 
nihilating light DM for DM mass fM = 0.75 MeV with 95% 
confidence constraint surfaces (dashed contours). Regions 
to the right of the constraint contours are ruled out. The 
rightmost constraint (orange) is ruled out by WMAP7+ACT 
01 and the middle constraint (green) by WMAP7+SPT @. 
The left constraint (blue) will be ruled out by forecast con- 
straints from Planck The dot indicates the Via Lactea 
II parameters. 



can also induce small distortions in the spectrum of the 
CMB (away from the Planck law) [tI-'so'I, which may be 
detectable in the future. Similarly, effects of energy de- 
position should be visible in the emerging 21 cm window 
for cosmological observations [sil - fssj . 



IV. 



CONSTRAINTS ON SCATTERING AND 
ANNIHILATING DARK MATTER 

A. Conservative Constraints 



The astute reader will notice a similarity between equa- 
tions © and (|9ll0p . which allows us, assuming that 
DM scattering is responsible for the e+ production in 
the galactic bulge, to find Pann or Pscatt as functions 
of the Einasto profile parameters a and r^. The key 
point to understand is the relation between {(Jann.iVrei) 
in Pann/ Pscatt, which is the thermal average in the early 
universe, to av. For s-wave annihilation, as appropriate 
for the MeV-mass annihilating DM of section III Al these 
quantities are equal. For XDM models with DM mass 
M > GeV, the DM velocity dispersion at recombination 
and later (but before structure formation) is considerably 
smaller than vt- Therefore, upscattering in endothermic 
XDM models is kinematically forbidden as noted, and 
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FIG. 3. Contours of logp^catt (in cm^/s/GeV) for XDM with 
vt = 10"^ and fSMi2 = 1.5 MeV with 95% confidence con- 
straint surfaces (dashed contours). Constraints are labeled as 
in figure (2] The dot indicates the Via Lactea II parameters. 



{criVrel) 



av for exothermic XDM since F{vrei = 0) = 1. 



As a result, we find that production of low-energy 
by exothermic scattering processes minimally deposits 



Pscatt — 



(11) 



This includes annihilation of light DM if we replace 
Ci — ^ Cann and set (5Mi2 = 2 A/0 As a result, constraints 
on Pann/Pscatt rulc out rcglons of a, Ts parameter space 
for light DM annihilation or exothermic XDM, assuming 
a fixed value of 5Mi2. Note that equation (fTTj) is indepen- 
dent of M except through a relatively weak dependence 
on Vt, which is determined by the ratio 5M2i/M . We will 
assume that / is constant in an on-the-spot approxima- 
tion; a more detailed parameter search using the redshift 
and energy dependence of / may be warranted in the 
future. 

Constraints on the allowed Einasto profile halo param- 
eters are shown in figure [2] for annihilating light dark 
matter and in figure |3] for exothermic XDM. For illustra- 
tive purposes, we have chosen 2fM,fSAIi2 = 1.5 MeV; 
this is a representative but conservative choice in that 
generated are mildly relativistic, lower e"*" energies 
require tuning in the DM model, and the constraints be- 
come more stringent for larger e"*" energies (M or dMi2). 
For reference, the Via Lactea II halo model for annihi- 
lating light DM is already ruled out by the 95% confi- 
dence constraints from WMAP7+SPT ^ if fM > 1.2 



MeV, and it is always ruled out by the forecast con- 
straints from Planck given in [t^ (unless the efficiency 
factor / is for some reason reduced below expected val- 
ues). As discussed earlier, we also choose vt — 10~^ in 
figure [3] as a conservative choice without fine tuning; for 
this Vt, the Via Lactea II model is already ruled out by 
WMAP7+SPT if SMi2 > 3.6 MeV (assuming / = 0.9) 
and would be free of the forecast Planck constraints if 
f5Mi2 < 1.2 MeV. For larger values of Vt, the constraints 
approach those shown in figure [21 

Similar constraints for annihilating light DM have been 
discussed previously in |26], which were based on the 
NF W profile selected by . While this profile is slightly 
cuspier in the innermost galaxy than the Via Lactea II 
Einasto profile, the corresponding value of Carm is ac- 
tually slightly less for the NEW profile assuming that 
both profiles are normalized to pq at tq. Furthermore, 
[l^ took p(7, = 0.3 g/cm^. Making the appropriate con- 
version, [26] constrains annihilating light DM to have 
/M < 1.7 MeV. As expected, the improvement in the 
CMB data since [26| has tightened the constraints. 

Finally, we note that these constraints are conserva- 
tive because they arise from the same mechanism that 
produces e"*" in the galactic bulge as long as the scat- 
tering is not velocity-suppressed (or kinetically forbidden 
at the time of recombination as in endothermic XDM). 
As a result, upcoming results from Planck will be able 
to exclude annihilating light DM and exothermic XDM 
at the Via Lactea II halo parameters (assuming that the 
mass splitting SM12 and threshold vt are not finely tuned 
in the case of XDM). In fact, the forecast Planck con- 
straints exclude a significant fraction of the best-fit pa- 
rameter space from [18j|. Aside from fine-tuning of the 
particle physics, our results suggest one additional way 
these models might evade CMB constraints. More cuspy 
Einasto profiles (with larger values of ^) will be uncon- 
strained by CMB experiments, and propagation of un- 
stable DM excited states or the produced e"^ themselves 
can, in principle, spread out the 511 keV gamma ray sig- 
nal. However, the similar morphology of the 511 keV 
and 130-135 GeV lines (pointed out in [i^) makes this 
way out less palatable, assuming that the 130 GeV signal 
withstands further scrutiny. 



B. Annihilation of XDM 

XDM can also annihilate completely, depositing addi- 
tional energy to the IGM, like a standard thermal WIMP. 
In the exothermic case, the total energy deposition is 
given by Pann + Pscatt in equation ([8]) . It is already pos- 
sible to rule out thermal WIMPs with mass in the range 
of 7 to 12 GeV, which has been of interest with respect 
to possible signals at the DAM A CoGeNT [ssl-lsej. 
and CRESST ,83| direct detection experiments, using the 
WMAP7+SPT constraints discussed above. In the fol- 



* Annihilation of DM produces one pair per DM pair vs one 
pair per excited DM state in XDM. 



lowing, we will illustrate that production of low-energy 
pairs provides extra sensitivity in the case of exother- 
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mic XDM in a model-dependent fashion. We focus on the 
exothermic XDM models of 0| as an example, working 
under a few simplifying assumptions. 

In these models, DM is a Majorana fermion triplet 
of a dark SU{2) gauge group and annihilates into the 
dark gauge bosons, assuming that the dark Higgs bosons 
are heavier than the DM. The dark gauge bosons all 
mix kinetically with the photon, so they can decay into 
any lighter charged SM particle. Then, for light gauge 
bosons, the final annihilation products are e^, /i^, and 
TT^. This allows us to find an average efficiency factor 
/ for the annihilation as the average of the "XDM elec- 
trons," "XDM muons," and "XDM pious" values of ^ 
weighted by the gauge boson branching ratios as a func- 
tion of the gauge boson mass. For simplicity and speci- 
ficity, we take all the gauge bosons to have the same mass 
of 500 MeV, which yields / = 0.53. 

In addition, the annihilation cross section at late times, 
{<^annVrei), differs from the canonical value 3 x 10~^® 
cm'^/s for three reasons. First, as discussed in [11], a more 
precise calculation finds a somewhat decreased value of 
the cross section needed for the correct thermal relic 
abundance in this mass range, and the required cross 
section has a significant dependence on mass. Second, 
as mentioned above, these XDM models contain dark 
gauge bosons which are mildly relativistic at the time 
the DM freezes out; as a result, there are more degrees of 
freedom in the primordial plasma than in standard cos- 
mology. This results in a slight decrease of the required 
cross section, as discussed in [13. Finally, we assume 
that the unstable DM excited state has decayed com- 
pletely by recombination (see section |V] below), so the 
average annihilation cross section changes after chemical 
freezeout because the DM states have different relative 
abundances. Specifically, at late times, the stable excited 
state has relative abundance Y while the ground state 
has relative abundance 1 — Y. Taking into account co- 
annihilations between the different DM states, {(JannVrei) 
is enhanced by a factor 98/75 at K = 1/3, its maximum 
value, in these SU{2) triplet models; this enhancement 
factor increases as Y decreases. We assume Y = 1/3. 

Results appear in figure |3] for DM mass M = 9 GeV 
and 10 GeV. We see that the exclusion region from the 
WMAP7-HACT constraint extends farther to the left (to 
cuspier halo parameters) than when not including DM 
annihilation (compare to figure [3]) . The improved sen- 
sitivity works in both ways: the WMAP7-I-ACT con- 
straint rules out these XDM models as thermal WIMPs 
for M < 8 GeV without XDM-like production, but 
the Via Lactea II halo is excluded for larger masses when 
both effects are taken into account. As a further exam- 
ple, the WMAP7-I-SPT constraint excludes these models 
for M < 18.4 GeV if XDM-like e± production is ignored 
but can exclude the Via Lactea II halo for M < 34.5 GeV 
when production is included. At the displayed DM 
masses, the annihilation cross sections at chemical freeze 
out are 2.39 x 10"^^ cm^/s for M = 9 GeV, 2.34 x lO'^^ 
cm^/s at M = 10 GeV, and 2.10 x IQ-^T cui^/s for M > 



15 GeV. 



V. CONSTRAINTS ON DECAYING DARK 
MATTER 

As discussed in section IIIIl the CMB also constrains 
energy deposition to the IGM from decay of DM. For r 
longer than the age of the universe, only a small fraction 
of the unstable DM has decays, so we find a constraint on 
Pdec- For lifetimes r shorter than the age of the universe, 
the exponentially decay of the abundance means that it 
is more convenient to constrain Yu{5Mi2/M) as a func- 
tion of T (and the decay product energy 6M12), where Yu 
now represents the initial abundance of the unstable DM 
species relative to all DM. |^ have provided constraints 
using WMAP3 data along with forecast constraints for 
Planck without estimating the deposition efficiency. Re- 
cently, |67|] updated these constraints to WMAP7 data 
(along with forecasting Planck constraints), including the 
dependence of the efficiency factor / on redshift and the 
decay product energy SM12. These constraints are in ba- 
sic agreement over a wide range of deposition energies, 
especially once a rough adjustment for the deposition ef- 
ficiency is included; we will use the results of [g^l, which 
are available as an electronic resource, as they use slightly 
more conservative assumptions. 

It is also important to note that the constraints are 
irrelevant for decays with lifetime much shorter than the 
time of recombination, ie t < 10^^ s, since faster de- 
cays deposit nearly all of their energy before the CMB 
decouples from the IGM. 



A. Decay of Unstable XDM State 

As noted in [s^l , XDM models necessarily have at least 
one unstable state (which we have denoted as state #2), 
and we can ask what constraints can be placed on the 
lifetime or initial abundance of that state. Here, we up- 
date the results of [33|. XDM may also have metastable 
excited states (denoted state #3), which must have a 
lifetime longer than the age of the universe, and we will 
consider constraints on this state, as well. 

We first consider the metastable state #3. For 
long lifetimes, the WMAP7 constraint becomes 
{SMi2/M){Y/t) < 10-24-8 s-\ while the forecast 
Planck constraint is {SMi2/M){Y /t) < 10^25.4 g-i^ 
Previous studies of XDM indicate that Y > 1/10 
with reasonable assumptions about kinetic freeze out, 
and it is furthermore difficult to arrange a sufficient 
scattering cross section for smaller relative abundance. 
Furthermore, it is reasonable to assume M > 10 GeV; 
with SM12 = 1.5 MeV, we find that r > 8.9 x 10^^ s 
(WMAP7) or T > 3.9 x 10^° s (Planck)- 

Turning to the unstable state #2, [i^l argued that the 
lifetime is constrained by the morphology of the 511 keV 
signal; if the lifetime is greater than approximately 10^"* 
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FIG. 4. Exclusion plots for a specific XDM model including energy deposition from DM annihilation and scattering for DM 
mass A/ = 9 GeV (left panel), 10 GeV (center), and 40 GeV (right). Contours are of pa„„ + Pscatt at the WMAP7+ACT 
constraint (orange) and the WMAP7+SPT constraint (green). Excluded regions are shaded. 
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FIG. 5. Maximum values of logcj = log(F„/Af) (in GeV~^) 
as a function of logr (in seconds) from WMAP7 (blue curve) 
and forecast for Planck (red dashed). This figure assumes 
5Mi2 = 1.5 MeV. 



s, the unstable DM will move outside the galactic cen- 
ter before emitting the pair. We show constraints 
from the CMB on oj = Yu/M for 5Mi2 = 1.5 MeV as a 
function of lifetime in figure [5l using the results of [i^l . 
The effect of the mass splitting 5Mi2 is mainly to shift 
the constraint curve vertically on the plot, but the shape 
of the curve also depends weakly on it. In most XDM 
models, Yu ^ Y due to the thermal history of DM; in 
this case, the lifetime can just be as long as 10^^ s if 
M > 19 GeV (WMAP7) or M > 77 GeV (Planck). For 
smaller masses, larger mass splittings, or larger initial 
relic abundances, we simply find r < 10^^ s. Other XDM 
models, such as in [43|, could have much smaller Yu and 
therefore potentially longer lifetimes. Note that these 
constraints are potentially stronger than the constraint 
from the morphology of the galactic signal but that the 
lifetime is degenerate with the mass, mass splitting, and 
relative abundance. These constraints apply in the case 
of either endothermic or exothermic XDM. 



We also considered the fact that a long lifetime for 
the unstable state could strengthen constraints on down- 
scattering in exothermic XDM models. Physically, if the 
unstable state is sufficiently long-lived, pairs are pro- 
duced at a later cosmological era than the initial down- 
scattering process. As a result, the DM density seems 
to be greater than expected in Pscatt- This process is 
governed by the Boltzmann equations 



dn2 
dt 



3Hn2 = 
3Hn3 ^ 



n2 



+ nl (aiVrel) 



-nliaiVrei) , 



(12) 



(ignoring kinetically forbidden upscattering) and energy 
deposition is governed by 



dE 
dVdt 



6Mi2n2{t) 



(13) 



The solutions to the Boltzmann equation for n2 can be 
given in terms of exponential integral functions (in the 
matter dominated era) , assuming the abundance of state 
#3 to be constant after early times. As expected, a life- 
time less than the time of recombination is equivalent to 
instantaneous decay. Meanwhile, a longer lifetime runs 
into the constraints discussed above unless the initial 
abundance of state #2 is quite small. In other words, 
this effect may be important in some region of parame- 
ter space for models like those of |43i]. 



B. Infalling Positron Model 



As discussed in section HTPl [20| recently proposed that 
decaying DM with lifetime in the range t = 10^^-10^1^ 
s and w ^ 4.7 x 10~^ GeV~^ produces sufficient e+ to 
explain the galactic 511 keV line. Note that w as defined 
in [13 is our % for M = 100 GeV; our w = Yu/M is 
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a measure of the number density of unstable DM par- 
ticlesH One point of importance is that the e"*" can be 
produced at high energies due to a mass sphtting of or- 
der GeV (or presumably more); they cool quickly in the 
IGM (mostly by Compton scattering) before entering the 
galaxy, where they thermalize and annihilate. But it is 
precisely the cooling process (and annihilation) in the 
IGM that modifies the CMB anisotropy spectrum. While 
detailed studies of the signal morphology for such mod- 
els have yet to be carried out, we show here that existing 
CMB constraints already rule out a great deal of the pa- 
rameter space of interest. 

Figure [S] summarizes constraints on these models. The 
blue band indicates values of u identified by [l^] as pro- 
ducing the correct number of e+ for the INTEGRAl sig- 
nal without overproducing them|f| The green region is 
excluded by limits on diffuse photons produced either by 
direct decays or Compton scattering of the resultant 
as found in [20| : this constraint lifts for mass splittings 
much below 1 GeV. The remaining regions are excluded 
by the WMAP7 constraints of [13 for 5Mi2 = 1 GeV, 
100 MeV, and 10 MeV (from bottom to top). We note 
that the preferred value of (5 A/12 = 1 GeV mentioned in 
[20| is nearly completely excluded: the central value of 
w = 4.7 X 10~^ GeV"^ lies above the WMAP7 constraint 
for T < 10^^ s and above the diffuse photon constraint 
for T > 10"'^'' s. Smaller values of the mass splitting 5Mi2 
open the parameter space somewhat, but part of the mo- 
tivation of the model is lost. Note that the CMB limits 
apply equally well to models in which the unstable DM 
state is charged and decays by emitting a single 6+ as 
opposed to an pair, since the e"*" can still efficiently 
deposit its kinetic energy and annihilate with an ambi- 
ent e~. As a result, we see that existing CMB data pro- 
vide robust constraints on these models. Larger mass 
splittings and the upcoming Planck data will make the 
constraints tighter. 

VI. SUMMARY 

As we have seen, CMB anisotropies provide robust con- 
straints on DM models for the galactic bulge e+ produc- 
tion mechanism. These constraints provide a new way to 
exclude annihilating DM and XDM models for a range 
of galactic DM halo parameters, including the fiducial 
parameters from the Via Lactea II simulation and much 
of the parameter space that best fits the 511 keV signal 
morphology [l^, as long at the scattering cross sections 
are not velocity-suppressed (or kinematically forbidden) 



^ Another difference is that we use the initial relative abundance 
of the unstable DM state, while [2Cll | uses the value at recombi- 
nation. Since the lifetimes considered are longer than the time 
of recombination, the difference is unimportant. 

^ We have estimated the width of this band, which is not specified 
in [131. 



-2 p 
-3 - 



-4 - 




log(r/s) 

FIG. 6. Constraints on infalling positron models in the 
logij /s) and log(tj GeV) plane. The blue band is the region 
of interest for the 511 keV signal. The top (green) region is 
excluded by constraints on diffuse photons for 5Mi2 > 1 GeV. 
Three lower regions from bottom to top are excluded by the 
WMAP7 constraints on decaying DM for 5Mi2 = 1 GeV (or- 
ange), 100 MeV (red), and 10 MeV (purple) respectively. 



at or after recombination and before structure formation. 
These models include s-wave annihilation as favored by 
the signal morphology [l^ and exothermic XDM models; 
while the Via Lactea II parameters are not yet excluded, 
Planck is expected either to exclude them or find signs 
of energy deposition. When including the effects of DM 
annihilation, the Via Lactea II parameters are in fact al- 
ready ruled out for some XDM models with DM masses 
less than a few tens of GeV, as we illustrated in one case. 

Our most conservative constraints arise because the 
same scattering events that generate the 511 keV gamma 
ray signal also produce e"*" in the early universe, which 
then deposit their energy in the IGM. The amount of 
energy deposition is set by the required rate of e"*" pro- 
duction in the galactic bulge. As a result, annihilating 
light DM and exothermic XDM cannot avoid these con- 
straints. Assuming the morphologically-preferred param- 
eter space is ruled out, these DM models will be viable 
only if the galactic DM halo profile is cuspier than usu- 
ally assumed. Then the morphology of the 511 keV signal 
would require that either the produced e"*" propagate sig- 
nificantly before annihilating or that the unstable excited 
DM state lives long enough to travel approximately 1 
kpc. Confirmation of the tentative 130-135 GeV gamma 
ray line at the galactic center would put pressure on this 
interpretation of the 511 keV signal, however, since the 
higher-energy gamma ray morphology appears consistent 
with the Via Lactea II halo. In this case, endothermic 
XDM could become the only remaining viable DM model 
for the 511 keV signal. A more optimistic point of view, 
on the other hand, is that Planck may provide evidence 
for nonstandard energy deposition in the IGM. In that 
case, more detailed studies will be warranted to deter- 
mine if that energy deposition is consistent with (or sug- 
gestive of) DM models for the galactic e+ production, as 
well as potentially measuring some parameters of those 
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models. 

In this paper, we also revisited constraints on the life- 
time of the excited DM states in XDM models by adapt- 
ing CMB limits on energy deposition from decaying DM. 
We further used these limits to exclude much of the pa- 
rameter space of interest in the recent proposal of [20| 
that 6+ produced by metastable decaying DM can gen- 
erate the 511 keV signal when they fall into the galactic 
bulge. 

In summary, the ability of DM models to explain the 
high rate of positron production in the galactic bulge, 
along with a striking and suggestive agreement between 
anticipated DM halo profiles and the signal morphology, 



motivates a search for additional signals of such DM. It 
is clear that CMB anisotropics provide a new handle on 
the viability of a broad class of these models along with 
their other parameters. 
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